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(54) Method and apparatus for identification of an access network by means of 1-port 
measurements 



(57) Method and apparatus for modeling a network 
are described in which a 1 port measurement (4) is 
made on the network by inputting an excitation signal at 
one port of the network and recording the results reflect- 
ed back to the port. A 1 port parametric model (6) of the — 
network is generated whereby as much Information 



about the topology of the network is included in calcu- 
lating initial values of the parameters. Then the values 
of the parameters of the 1 port parametric model are 
optimized (7) by reducing the difference between the re- 
sults of the measurement step and resu Its calculated us- 
ing the 1 -port parametric model and the excitation sig- 
nal. 
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Description 

[0001] The present invention relates to simulation methods for telecommunications networks and to telecommuni- 
cations networks themselves having a receiver and a transmitter connected by some form of communications channel. 
5 The present invention is especially useful for the prediction of network parameters such as a capacity of the network, 
and in particular xDSL networks. 

Technical Background 

^0 [0002] The acronym xDSL stands for the family of Digital Subscriber Line technologies, which allow high-speed 

access to the Internet and multimedia services over the local loop, which connects the CP (customer premises) to the 
CO (central office), that is over simple twisted pair cables. An xDSL transceiver at the CO communicates with an xDSL 
transceiver at the CP over the local loop. 

[0003] Since decades the local loop, which is a transmission line consisting of two twisted copper wires, also called 
unshielded twisted pair (LTTP), has given the customer access to POTS (Plain Old Telephony Service). The POTS 
signal, transmitted over the local loop, is analog and contained in the frequency band up to 4 kHz, which corresponds 
to the spectral content of speech. 

[0004] xDSL exploits the frequency band above 4 kHz up to several MHz, which is not used by POTS. However as 
the legacy local loops have been engineered for voice-band transmission, there are no guarantees about the quality 
^ of ttic local loop with respect to transmission in this higher frequency band. Consequently it is preferable If every loop 
IS qualified for xDSL. 

[0005] The quality of a loop is expressed by Its theoretical channel capacity, which is equivalent to the upper bound 
of the achievable bit rate. The capacity depends on the signal-to-noise ratio (SNR) as a function of the frequency at 
the receiver at the CP, respectively the CO, for the downstream transmission (from the CO to the CP), respectively the 
upslroam transmission (from the CP to the CO). The SNR at the receiver at the CP. respectively the CO, is ‘determined 
as a function of the frequency by the transfer function of the loop between the CO and the CP and the noise PSD 
(Power Spectral Density)at the CP, respectively the CO given the PSD of the transmitted signal at the CO, respectively 
the CP 
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[0006] In general the local loop consists of a network of transmission lines. Every fine in the network is a LTTP char- 
acterized by its length and type. The line type specifies the cross-sectional geometrical dimensions, such as the wire 
diameter (also called wire gauge), and the material physical constants, such as the electrical permittivity of the dielectric 
separating the 2 copper wires. The most used wire diameters are 0.4 mm, 0.5 mm, and 0.6 mm. Polyethylene (PE) is 
the most occurring insulator but other materials are also used such as paper and PVC. 

[0007] The network topology of the local loop is limited to a tree structure. The simplest topology is a single line. The 
magnitude of the transfer function reflects the attenuation of the line, which gets worse with increasing frequency and 
line length. Another topology that exists for long loops is a cascade of 2 or more lines with increasing wire diameter 
from the CO to the CP. For this topology reflections are caused by the change of the wire diameter at the splices 
connecting 2 lines. A topology that is also frequently encountered, especially in the USA, is a loop with 1, 2 or more 
bridged taps. A bridged tap is an openended line spliced to the main line. Reflections appear for this topology at the 
splice connecting the bridged tap to the loop and at the open end of the bridged tap. Reflections have a negative Impact 
on the transfer function, because they interfere with the signal propagating along the direct path. For those frequencies 
for which the interference is destructive, the magnitude of the transfer function reduces. Such reductions rarely appear 
in the voice band because the bridged taps are usually not too long. 

[0008] To improve the voice-band transmission for very long loops it has been a common practice to install loading 
coils, which are inductors of typically 88 mH, inserted in series in the loop. Typically the first loading coil is at 900 m 
from the CO and subsequent loading coils are spaced 1800 m apart. These loading coils act as low-pass filters with 
4 kHz as cut-off frequency, which has such a bad effect on the transfer function above 4 kHz that xDSL cannot operate. 
[0009] As the twisted pairs constituting the local loop are unshielded, external electromagnetic waves may couple 
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into the loop and propagate towards the CO and the CP causing noise at the receiver. The electromagnetic coupling 
is reduced by the twisting of the 2 wires, because adjacent segments of the twisted pair experience electromagnetic 
waves with opposite polarity. In addition the twisting improves the balance of the line. A line is balanced when the 2 
conductors have an equal impedance towards the earth. The balancing of the line prevents a common-mode signal 
from transforming into a differential-mode signal. In the case of a common-mode signal the 2 wires carry equal currents 
and the return path of the current is the ground. For a differential-mode signal the 2 wires carry opposite currents (out 
of phase currents). Electromagnetic waves may couple into the line because of imperfect twisting, and the common- 
mode signal that they cause, may transfonm into a differential-mode signal because of imperfect balancing, which is 
correlated with the twisting. Balance decreases with increasing frequency. 

[0010] The noise is divided into 2 different types according to the origin of the external electromagnetic waves cou- 
pling into the loop. The first type of noise is crosstalk, which is the electromagnetic coupling between twisted pairs in 
the same cable. The cables leaving the CO contain thousands of twisted pairs. The closer to the CP the less pairs 
there are present in a cable. A difference is made between near-end crosstalk (NEXT) and far-end crosstalk (FEXT). 
The transmitters at the CO. respectively CP, are the source of NEXT for the near-end receivers at the CO, respectively 
the CP, and are the source of FEXT for the far-end receivers at the CP, respectively the CO. In general crosstalk gets 
worse with increasing frequency. 

[0011] A second type of noise is radio-frequency interference (RFI), which is caused by radio waves coupling Into 
the local loop, that acts as an antenna, especially if there are aerial lines. There are 2 major sources of radio waves 
in the frequency band of xDSL namely AM radio and amateur radio. 

[0012] Hence, the local loop has several impaimients for transmission In the frequency band of xDSL, which are not 
present for voice-band transmission. The existing metallic line testing (MLT) systems for POTS are not capable of 
qualifying the loops for xDSL accurately. Telecom operators need a reliable qualification tool for the deployment of xDSL. 
[0013] For example, if a customer requests a certain xDSL service, the operator has to be able to estimate the 
achievable bit rate for the loop that connects the CP of that customer to the CO, especially when a precisely specified 
quality of service (QoS) in terms of bit rate is to be guaranteed. If the achievable bit rate is underestimated, there is 
the risk of lost business. If the bit rate is overestimated, the telecom operator risks ending up with a dissatisfied customer 
and a useless Installation as well as troubleshooting costs. As the latter risks are the most Important, the telecom 
operator will subtract a safety margin from the estimated bit rate. In this way the risk of loosing customers and having 
non-refunded costs is reduced at the expense of increasing the risk of lost business. The more accurate the estimate 
of the bit rate is , the smaller the safety margin can be, and the less lost business. Hence there is a need for an estimation 
of the bit rate which should be as accurate as possible, but without incurring large additional costs. 

[0014] Currently there are systems available on the market for xDSL loop testing, but they require measurements 
at the CO and at the CP. So a technician has to visit the CP as well as the CO. Furthermore, the interpretation of the 
measurement results, which requires a lot of expertise, is mainly left to the technician. This makes the loop qualification 
labour intensive and hence too expensive, in general. Telecom operators are not eager to use this type of qualification 
scheme on a large scale. The way the loop qualification Is actually handled d^ends on the operator. If there is a 
database available containing information about the loops in the access network, this knowledge can be used for the 
loop qualification. If the topology and the lengths and types of the individual line sections of the loops are available, 
then the transfer functions of the loops can be computed. In order to determine the theoretical channel capacity, the 
noise PSD at the CO and at the CP is still missing, but the transfer function is already very infomriative. However not 
all telecom operators maintain such a database. Furthermore, if the database exists, It often contains a lot of errors, 
because it is difficult to keep it up-to-date. Without a database a telecom operator has to make do with rules of thumb. 
For example the distance from the CP to the CO can be used as very rough criterium for loop qualification. The con- 
clusion is that telecom operators need a completely automated xDSL loop qualification system that performs meas- 
urements at the CO end of a loop and can interpret the results in order to estimate the achievable bit rate of the loop 
as accurately as possible. When the loop qualification has been finalized, this system can be used as an xDSL loop 
testing system for maintenance purposes. 

[0015] In order to detemiine the achievable bit rate of a loop, the transfer function of the loop and the noise PSD at 
the CP, respectively the CO, are needed for the downstream, respectively upstream transmission. These quantities 
have to be measured from the CO end of the loop. For the noise PSD at the CO this is trivial. The loop transfer function 
and the nose PSD at the CP cannot be measured directly at the CO. 

[0016] The quantity that is directly measurable at the CO, which is named port 1 , is the scattering parameter S11 , 
which is the ratio of the reflected wave and the incident wave at port 1 . The reflected wave is also called a reflectogram 
and the measurement procedure is named reflectometry, more specifically time-domain reflectometry (TDR) if the 
waves are measured in the time domain. As the loop Is considered as a 1-port, the scattering parameter S11 quals 
the reflection factor and the loop is completely characterized. However as the loop is designed for transmission between 
a CO and a CP, It has to be considered as 2-port with the CO as port 1 and the CP as port 2. Then the loop is completely 
identified by its 4 scattering parameters S11 , S21 , S12, and S22. The transfer function is related to the S21 scattering 
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parameter. Time-domain reflectometry (TDR) is a well known measurement technique, which has been applied since 
a.few decades in different domains such as cable fault location. Since the beginning of the xDSL deployment TDR has 
been regarded as a candidate technology to solve the loop qualification problem. There are now xDSL loop testing 
systems available with an integrated TDR meter, but the interpretation of the measured reflectogram is left to the 
technician. The loop qualification performed with these systems aims at the detection, location, and removal of local 
loop impairments such as bridged taps and load coils. 

[0017] Recently attempts have been made to process the reflectogram in order to estimate the theoretical channel 
capacity of the loop. This processing is done by means of a artificial neural network (ANN), which transforms a number 
of inputs into a number of outputs by means of some elementary mathematical operations such as addition or multi- 
plication with scalar, and/or non-linear functions. These operations are structured in a way that resembles a real neural 
network. The specific structure and the values of the multiplying scalars, called weighting factors, determine the input/ 
output behaviour of the ANN. 

[001 8] Some features of the reflectogram such as the position and height of the peaks, corresponding to the reflected 
pulses, are used as inputs and some interesting parameters of the loop are used as outputs such as the attenuation 
of the loop at a number of frequencies. 

[0019] There are a number of disadvantages related to the use of an ANN network. The ANN has to be trained, I.e. 
its structure has to be defined and the weighting factors have to be tuned, using a large number of sets of inputs and 
corresponding outputs. The training procedure has to cover all combinations of loop and noise types present in the 
network. Extrapolation with respect to the training set is unreliable. For example If the training has been done for loops 
up to 2 km, then the neural network will give an unpredictable output for a loop of 3 km. Only interpolation Is allowed, 
but validation is necessary. Due to the increasing complexity of the ANN with an increasing number of input and output 
arguments, the number of inputs and outputs of the ANN has to be limited. This means that not all the information 
present in the reflectogram is used and this has a negative impact on the accuracy. In addition only a few parameters 
related to the loop can be predicted. An ANN is a black box model. It doesn't allow Inclusion of some priori knowledge 
about the loop that is available, e.g. the fact that a loop consists of a network of transmission lines. This lowers the 
accuracy. 

Summary of the Invention 

[0020] An aspect of the present invention may be descnbed as a method and apparatus for the determination of a 
model of a 2-port system based on a 1 -port measurement and an optimization of a 1-port parametric model using the 
measurement. The optimized parameters are used in a 2-port parametric model of the system in order to determine 
useful characteristics of the system such as transmission capacity between the two ports. 

[0021] The present invention also provides a method and apparatus for determining the expected capacity of a 
transmission line only using the measurement of the scattering parameter for which the transmission line is con- 
sidered as a 1-port, andthe measurement of the noise PSD at one end of the transmission line, e.g. in a GO -Capacity 
estimation is only one possible application of the present invention. 

[0022] Another aspect of the present invention is the estimation of the end-to-end transfer function of at least one 
transmission line based on measurements at one termination location thereof only. It Is not possible to derive the 
transfer function directly from such single-ended measurements. However measurements at one end of the loop allow 
characterization of the loop completely. The result of the characterization is a model of the transmission line, which 
can be used to predict the transfer function. 

[0023] In order to characterize the transmission line an excitation signal is generated at one termination location 
thereof, e.g. at a central-office end of a local loop, and at the same location the response signal is measured. The 
excitation signal corresponds to the incident wave and the response signal to the reflected wave. The latter is also 
called reflectogram. This measurement procedure is called reflectometry (time-domain reflectometry, or TDR, if the 
signals are measured in the time domain). The ratio of the reflected wave and the incident wave is defined as the 
scattering parameter S^., (also called reflection factor in the case of a 1-port). So the scattering parameter of a 
certain transmission line describes how that line responds in terms of a reflected wave, if it is excited with a certain 
incident wave. It is remarkable that it is possible to model such a transmission line completely, including its transfer 
function, only measuring the scattering parameter complete model of the line considered as a 2-port consists 

of the 4 scattering parameters S 21 , S| 2 i snd The transfer function is related to the scattering parameter S 21 . 
[0024] In one aspect the present invention makes use of a 1 -port parametric model of the at least one transmission 
line in addition to the scattering parameter measurement to det rmine a 2-port parametric model of the at least 
one transmission line. A transmission loop, e.g. a local loop usually consists of a network of transmission lines con- 
necting the customer premises to the central office. The network topology is limited to a tree structure. Every line in 
the network is characterized by its length and type. The line type specifi s the cross-sectional geometrical dimensions 
and the material constants. 
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[0025] The present invention may provide a method to model a representation of a network of at least one transmis- 
sion line, said method comprising the steps of: 



generating and sending over said at least one transmission line an excitation signal from one end of said at least 
5 one transmission line; and 

receiving at one end of said at least one transmission line a reflectogram and post-processing said reflectogram 
to form a frequency domain reflectogram; 

wherein said post-processing includes the steps of: 
w 

generating a modeled frequency domain reflectogram using a frequency domain parametric model of said at least 
one transmission line having adaptable parameters; 

quantifying a difference between said modeled frequency domain reflectogram and said measured frequency do- 
main reflectogram; and 

- reducing the difference by adapting parameter values of said modelled frequency domain reflectogram. 



[0026] The method may also comprise the step of: 

Fourier transforming said measured reflectogram into said measured frequency domain reflectogram; 

[0027] The present invention may also provide a method to model a representation of at least one transmission line, 
said method comprising the steps of: 

generating and sending over said at least one transmission line an excitation signal from one end of said at least 
one transmission line; and 

receiving at one end of said at least one transmission line a time domain reflectogram and post-processing the 
reflectogram; 



wherein said post-processing includes the steps of: 

30 

generating a modeled time domain reflectogram using a frequency domain parametric model of said at least one 
transmission line having adaptable parameters; 

quantifying a difference between said modeled time domain reflectogram and said measured time domain reflec- 
togram; and 

35 - reducing the difference by adapting parameter values of sard frequency domain parametric model. 

[0028] In the above methods the post-processing step may include estimating a transfer function of said at least one 
transmission line. 

[0029] The step of generating a modeled time domain reflectogram In the above methods may include: 

40 

generating the frequency domain parametric model for said at least one transmission line; 
transforming said frequency domain parametric model into the modeled frequency domain reflectogram. 

[0030] The above methods may include a step of reducing the difference which Includes: 

45 

quantifying said difference as a cost function; and 



iteratively minimizing said cost function by adapting parameter values of said frequency domain parametric model. 
The cost function may be a maximum likelihood estimator or a Bayesian estimator. 

[0031] The output of any of the above methods can be at least one of: 



a set of optimized parameters for the model which minimizes the difference 
a transfer function of the at least one transmission line 
an estimated reflectogram based on the optimized parametric model 
a comparison between an optimized reflectogram and the measured reflectogram 
an estimated capacity for the at least one transmission line 

a yes/no answer or a probability that the at least one transmission line can support a capacity abov or below a 
predetermined target capacity 
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a message that the at least one transmission line cannot be modeled within a predetermined accuracy 
-- a length of the at least one transmission line which can support a predetemnined capacity. 

[0032] In the above methods the at least one transmission line can be a local loop of an access network. 

[0033] The present invention may also provide a system for modelling a representation of at least one transmission 
line, comprising: 

means for generating and sending over said at least one transmission line from one end thereof an excitation 
signal; and 

means for receiving at one end of said at least one transmission line a reflectogram and means for post-processing 
the reflectogram to a frequency domain reflectogram, 

wherein said means for post-processing said reflectogram comprises: 

means for generating a modelled frequency domain reflectogram using a frequency domain parametric model of 
said at least one transmission line having adaptable parameters; 

means for quantifying a difference between said modelled frequency domain reflectogram and said measured 
frequency domain reflectogram; and 

means for reduang the difference by adapting parameter values of said modelled frequency domain reflectogram. 

[0034] The present invention may also provide a system for modelling a representation of at least one transmission 
line, comprising: 



means generating and sending over said at least one transmission line an excitation signal from one end of said 
at least one transmission line; and 

receiving at one end of said at least one transmission line a time domain reflectogram and post-processing the 
reflectogram; 



wherein said means for post-processing includes: 

means for generating a modelled time domain reflectogram using a frequency domain parametric model of said 
at least one transmission line having adaptable parameters; 

- means for quantifying a difference between said modelled time domain reflectogram and said measured time 
domain reflectogram; and 

means for reducing the difference by adapting parameter values of said frequency domain parametric model. 
[0035] The means for generating a modelled frequency domain reflectogram may comprise: 



means for generating the frequency domain parametric model for said at least one transmission line; 

means for transforming said frequency domain parametric model into a modelled frequency domain reflecto- 
gram; and 

Fourier transforming means for transforming said measured time domain reflectogram into a measured fre- 
quency domain reflectogram. 

[0036] In the above systems the processing means may include means for estimating a transfer function of said at 
least one transmission line. 

[0037] In the above systems, the means for reducing the difference may include: 
means for quantifying said difference as a cost function; and 

means for iteratively minimizing said cost function by adapting parameter values of said frequency domain para- 
metric model. 



[0038] The above systems may include means for outputting at least one of: 

55 

- a set of optimized parameters for the model which minimizes th difference 
a transfer function of the at least one transmission line 
an estimated reflectogram based on the optimized parametric model 
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a comparison between an optimized reflectogram and the measured reflectogram 
- - an estimated capacity for the at least one transmission line 

a yes/no answer or a probability that the at least one transmission line can support a capacity above or below a 
predetermined target capacity 

5 - a message that the at least one transmission line cannot be modeled within a predetermined accuracy 

a length of the at least one transmission line which can support a predetermined capacity. 

In the above systems the at least one transmission line may be a local loop of an access network. 

[0039J The above systems for modeling may be included in a central office system of a telecommunications network. 

10 [0040] The present invention also includes a data carrier medium carrying one or more computer readable code 

segments for controlling a processing device to carry out any method in accordance with the present invention. 
[0041] The present invention also includes a computer program product for execution of any of the methods in ac- 
cordance with the present invention on a computer system. 

[0042] The present invention also includes a method for modelling a representation of network, said method com- 
prising: 



transmitting a description of the network from a near location to a remote computing system via a telecommuni- 
cations network, 

executing on the remote computer system any of the methods of the present invention, and 
transmitting to a near location information relating to the response of the network to an excitation signal. 

The present invention will now be described with reference to the following drawings. 

Brief Description of the Drawings 
[0043] 



Fig. 1 trartsmission line model of a single loop. 

Fig. 2 transmission line model of a bridge tapped loop. 

Fig. 3 transmission line model of a loop of two cascaded lines. 

Fig. 4 measurements and estimated values of the scattering parameter versus frequency using embodiments 
of the present invention for a single line. 

Fig. 5 measurements and estimated values of the transfer function H versus frequency using embodiments of the 
present invention for a single line. 

Fig. 6 measurements and estimated values of the scattering parameter S.|^ versus frequency using embodiments 

of the present invention for- two cascaded lines. * • 

Fig. 7 measurements and estimated values of the transfer function H versus frequency using embodiments of the 
present invention for two cascaded lines. 

Fig. 8 measurements and estimated values of the scattering parameter versus frequency using embodiments 
of the present invention for a bridge tapped line. 

Fig. 9 measurements and estimated values of the transfer function H versus frequency using embodiments of the 
present invention for a bridge tapped line. 

Fig. 10 is a schematic representation of a computing device which can be used with the present invention. 

Fig. 11 is a schematic representation of a method flow diagram of an embodiment of the present invention. 

Description of the details embodiments 

[0044] The present invention will be described with reference to certain embodiments and drawings but the present 
invention is not limited thereto but only by the claims. In the following emfc>odiments of the present invention will be 
described by a parametric model detailed for a limited number of network topologies but this does not limit the appli- 
cability of the invention. The skilled person will appreciate that the present invention may be extended to other topologies 
and to other materials once the principles of the present invention have been understood. All such additional topologies 
are included within the scope of the present invention. Also in the description a twisted pair transmission line model is 
given. However the invention is not limited to this model and is not limited to twisted pair cables as transmission lines. 
The application of the pres nt invention is not limited in this respect and it can be used for symmetric or asymmetric 
lines, for example in the cas of coaxial cables, optical fibers, etc. and even in the cas of a general 2-poit network. 
[0045] The following references may be found useful for understanding the present invention. 
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[1] P. Boets, "Frequency Domain Identification of Transmission Lines from Time Domain Measurements", Ph-D. 
Thesis, Vrije Universiteit Brussel, Dept. ELEC, Pleinlaan 2, 1050 Brussels. June 1997. This document is Incorpo- 
rated by reference in its totality. 

[2] P. Boets, M. Zekri, L. Van Biesen, T. Bostoen, andT. Pollet, "On the Identification of Cables for Metallic Access 
Networks," in Proc. IMTC, 2001, incorporated by reference in its totality. 

[3] D. M. Pozar, Microwave Engineering, 2nd ed. New York: John Wiley & Sons, 1998. T Starr, J. M. Cioffi, and P. 
J. Silverman, Understanding Digital Subscriber Line Technology 
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[0046] In accordance with an embodiment of the present invention, given a specified or assumed network topology, 
a 1 -port parametric model for the scattering parameter of the network is derived from the models for individual line 
sections of the network or based on an assumed model. In an S.,., scattering parameter model the parameters of the 
individual lines appear as independently identifiable parameters. It is then possible to model the loop completely by 
only measuring the scattering parameter S^^ . The unknown parameters have to be estimated by a further optimization 
routine. The optimized parameters are used in a 2-port model of the network for further analysis of the network. 
[0047] Accordingly, each line in a network is fully characterized by a parametric model, which consists of the char- 
acteristic impedance and the transfer function of the line as a function of the frequency and a number of unknown 
parameters. Several line models are available. A model based on electromagnetic theory is one of the most useful for 
the line characterization. This model is based on the 4 primary transmission line parameters, which are resistance, 
inductance, conductance, and capacitance, as a function of the frequency and a number of parameters related to the 
cross-sectional geometrical dimensions and material constants. These latter parameters are considered to be un- 
known. In the following, it will be assumed that the transmission line comprises twisted pair cables. For other cable 
types, equivalent models may be selected. For twisted pair cable, the electromagnetic model assumes that the wire 
for the local loop comprises 2 cylindrical conductors with the same diameter immersed in a homogeneous dielectric. 
It allows for the skin effect and the proximity effect (up to 2 terms in the Fourier series solution). It contains 5 independ- 
ently identifiable parameters, which are related to the cross-sectional geometric dimensions and material constants of 
the components of the twisted pair cables. The electrical permittivity of the dielectric is considered to be independent 
of the frequency in this model and the conductance is neglected, which are good approximations for polyethylene (PE). 
This line model will be called VUB1 . 

[0048] Parameter estimation is an important part of characterization of a loop in accordance with the present inven- 
tion. In accordance with an embodiment of the present invention the difference between measured values for a loop 
and the modeled values for the same loop is reduced by appropriate selection/deteimination of the unknown param- 
eters. For example, a cost funetbn may be introduced which quantifies the deviation between the measured values 
(reflectogram) and the result of modeling the scattering parameter of the loop as a function of the model parameters. 

In order to match the model as well as possible to the measured values the cost function Is minimized as a function of 
the model parameters. The values of the model parameters corresponding to this minimum are the estimated values 
of the parameters. The lower limit of the variance of the estimated parameters may given by any surtable.method, e. 
g. by the Cramer-Rao lower bound. The cost function is also called an estimator, because it determines the parameter 
stimation completely. Several estimators are available. The maximum likelihood estimator (MLE) is preferred because 
it is consistent, asymptotically unbiased, asymptotically normally distributed, and asymptotically efficient under certain 
noise assumptions. In addition it is robust. It has such good properties because it takes the noise in the measurement 
into account. This means of course that this noise has to be measured. The characterization is typically performed in 
the frequency domain and the measurement in the time domain but the present invention is not limited thereto. Meas- 
urements and characterization may be performed in either the time domain or the frequency domain. 

[0049] The minimum of the cost function as a function of the model parameters is found by means of an optimization 
algorithm. Typically, this optimization algorithm will be iterative There are several known optimization routines which 
may be used. The Levenberg-Marquardt (LM) method Is preferred. It combines the Gauss-Newton and gradient descent 
procedures. The Gauss-Newton method is very fast, because it makes use of second-order derivatives, but it does 
not always converge. This problem of convergence is solved in the Levenberg-Marquardt algorithm by addition of the 
gradient descent method. The Levenberg-Marquardt algorithm only searches for a local minimum. In order to find the 
global minimimum, the initial values of the model parameters are preferably chosen to be close enough to the exact 
values, which correspond to the global minimum of the cost function so that the optimization does not become trapped 
in a local mimimum. As the estimated parameters are available for all the individual lines, the loop is completely modeled 
and the estimated transfer function can be calculated. 

[0050] In accordance with an embodiment of the present invention, a suitable optimization is based on aS.|^ scattering 
parameter measurement, an MLE cost function and an LM optimization algorithm. The scattering parameter model 

depends on the network topology and the model for the individual line sections. For every different network topology 
and line model there is a different S.,., scattering parameter model. This model has to be calculated. Independently 
identifiable parameters have to be determined. The derivatives of the model to its parameters have to be computed. 



8 



OCIO: <EP 1248363A1 I > 




EP 1 248 383 A1 



because they are needed for the Jacobian in the LM optimization algorithm. 

[0051] This has been illustrated for 3 different loop topologies: a single line, a cascade of 2 lines, and a loop with 1 
derivation, also called bridged tap. but the skilled person will appreciate after having understood the present invention 
that i: may be extended to other topologies. For all these loop topologies the assumption has been made that the loop 
5 IS open at the customer-premises end, which is a fair approximation of a phone set that Is on-hook. For the loop with 
a Dridgod tap an additional assumption is that all the individual lines are of the same type. These 3 loop topologies 
hrtve been chosen to cover all the different causes for reflections In a loop: the mismatch between the generator 
irrpodance and the characteristic impedance of the line connected to it, the mismatch between the load impedance 
and nc characteristic impedance of the line connected to it, the mismatch between 2 lines with a different characteristic 
to irrpcdancc connected to each other, and the mismatch due to the connection of a derivation. 

[0052] When the topology of a loop Is unknown an alternative method in accordance with an embodiment of the 
present invention tries each one of a series of well known topologies and makes an optimization of the unknown 
pr-.rameters lo' each topology by comparison with the reflectogram. The topology which appears to give the best fit 
with the ref cc:og-am may then taken as the best fit topology for the calculation of the network capacity. 

[0053] A base methodology in accordance with the present invention will be described with reference to Fig. 11. 
^ rsii^ a 1 port paramotnc model 2 having a finite number of coefficients whose values need to be optimised is generated 
m step ' fer a network based on an actual or assumed topolgy 3 for the network. Then a 1-port measurement 4 is 
earned out in step 5 on the network to be analysed. For example, a measurement of the scattering parameter may 
be carried out but the present invention is not limited thereto - the input impedance Zin could equally well be used. 
• Ttiib nieaj>cren»crit is used to do the parameter optimisation of the 1 -port model 6 of the network in step 7. Then these 
cstrraied parameters 9 can be included within a 2-port parametric model 1 0 for the access network based on a trans- 
m ssion line model lor every Imc segment to compute the estimated transfer function H in step 8. For the purposes of 
v«tiKtitioo only this end to end transfer function can be validated with the measured transfer function and the transfer 
funct on cs: mated direct ty via ihc transfer function measurement, however, this is only for valmidation purposes. 
(00S4J The 2 port model of Iho network can then be used for analysis pruposes, e.g. for the calculation of capacity. 
[OOSS] T ho prosoni invont*on is not limited to access networks of telephone networks, it can be used for any medium 
used for If Mhc'c there is a transmitter - medium - receiver structure. To identify the medium, 1 port meas- 

uromonis h! tr ^nvn^tor or receiver side can be made and these used to optimisze a parametric 1-port model for the 
modi jm T he .metre frxKdci preferably has a limited set of parameters. The parameters of the transmission medium 
^ mDdci arc opt m,voc ard then used to determine a 2-port model for the network. This 2-port model can be used for 
aralysis of tra-*smisk*on characteristics. The transmission medium could be any one suitable for a coimmunications 
system c g a r>o?wor% of twisted pair lines, coaxial lines, optical fibers, the air interface of a mobile teleocmmunications 
network 

35 I MODELS 



[0056] In acco^darKTc with embodiments of the present invention a transmission line or a combination or network of 
transmiss»cn itocv a rrx>dc«od by a parametric model having a finite number of coefficients to be determined, e.g. 5 
coeff cionis to lor each "oasic transmission line" or transmission line primitive. Such a transmission line primitive 
can conrprisc a rwtsted pair cable such as is used conventionally for the local loop of a public telephone system or, for 
instance a iw sled pair connection between two network devices on a Local Area Network. To model a network of 
transmission ir^cs each transmission line primitive is represented by the above parametric model - that is the network 
IS mcdo'od by a set of individual parametric models, one of each transmission line primitive and each having a finite 
number of coctcicnis to be determined. Some of the coefficients (e.g. to a^) may be estimated initially by calculating 
an in iiai vaiiio based on characteristics of the relevant transmission line, e.g. physical constants of materials used in 
the conrcctmc; c^»c forminq ir>o transmission line. Other coefficients may be selected, e.g. a^., however it is preferred 
if ihc:>o coeff«ncrKs arc also oeicf mined SO as to bring the initial model close to the optimized one. Knowledge of the 
actual topo o' trx? ixrtworK or of a finite number of topologies which the network may have can be used to set up 
the parametric model 

[0057] In H f jrihcf step tic reflectogram of the network is measured experimentally from one access point or "port" 
of the ncivkrorK To rcco'd the reflectogram a suitable test wave or pulse is entered at the relevant port and reflected 
waves mc.isufcd at :nc port. This reflectogram is preferably stored in a manner that allows manipulation of the reflec- 
togram within ^ a qiiai computing device, e.g. the reflectogram may be converted to digital signals using an analogue 
to digitai converter In a next step the finite number of coefficients of the paramteric model are optimized so as to give 
a close (or at least a c oser lit than the initial values) between a characteristic of the network as modeled and the 
reflectogram This rompar son may be done in the time domain or the frequency domain. In the following only a com- 
parison in the Ifcqiicrvry domain wii: be described. For example, an optimization routine may be used to optimize the 
coefficients of the paramclrc model The optimization may be done with respect to a cost function. That is, the goodness 
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of fit between the reflectogram and the results predicted by the parametric model (e g. a difference) may be represented 
as a quantitative value, and the coefficients may be optimized with reference to this value. For instance, this value may 
be maximized or minimized depending on whether the chosen value is a large or small value, respectively when there 
is a good fit. 

[0058] One problem with such an optimization is that the routine may stop in a local minimum or maximum and not 
an absolute one. To avoid or reduce the frequency of this problem the present invention makes use of a quasi -realistic 
model of the network so that there is some correlation between at least certain coefficients of the parametric model 
and physical constants of the network. By this means at least some of initial values of the parameters should be close 
to their final optimized ones and therefore there is less chance that the optimization routine lands in a local minimum. 
This is a significant advantage of the present invention compared with fully heuristic or pragmatic models for which 
initial values cannot be calculated reliably from network materials and topologies. 

[0059] Finally, the optimized parameters may be used in a 2-port parametric model of the network to determine 
characteristics thereof, e.g. a capacity of the network, more specifically a capacity which provides a certain Quality of 
Service (QoS) such as a threshold signal to noise ratio (SNR), Bit Error Rate (BER), Symbol Error Rate (SER), Frame 
Error Rate (FER) or similar as is known to the skilled person. A noise power spectral density (PSD) at the ends of the 
network, e.g. the CP and CO, is required for these calculations (see formula based on the theory of Shannon given as 
equation 1 above for capacity determination in bits per second). It is generally always possible to measure the noise 
at one location or port. The noise at the other location or port has to be estimated or a model has to be used. 

[0060] In the following, a non-limiting number of models of useful topologies for local loops of telephone systems will 
be described. 

A. Transfer Function H and Characteristic impedance 2J, ot an Unshieided Twisted Pair (UTP) Transmission Line 

[0061] The following equations may be used for describing a basic parametric model in accordance with an embod- 
iment of the present invention for an unshielded twisted pair transmission line. This basic parametric model will be 
used to construct more complex models of networks which can be described by a combination of basic parametric 
models. 



which represents the propagation constant of the line times its length. 

* • ■ - 3 — - — — ~ 2 - • 

a^d^ «y2+3a2*^o*/3+'^o*^i 

which defines an auxiliary function used in the calculations. 



with /= 0,1, 2, 3 (J = Bessel function) 






ry^ 



(5) 

( 6 ) 



which represents the transfer function of a primitive transmission line section 



yp/= sgs 

where yp represents the parallel per-unit length admittance of the transmission line 




(7) 



( 8 ) 
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where Zq represents the characteristic impedance of the line 

[0062] In this model there are five coefficients which need to be optimized with respect to a measured reflectogram 
The coefficients to be determined (a, to a^) may be given initial values in accordance with: 
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• air Al rr 



an^io 



(9) 
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( 10 ) 



a^ = aj^ ( 11 ) 

15 

a 

^ ln(1/7i^) 

i- The initial value of is prefereibly chosen as a combination of initial values of a-,, 82 . and 83 as explained in ref. 1 
which is incorporated herewith by reference. 



as=C/ (12) 

25 

[0063] By introducing parameter a^^ no reliance should be placed on the physical meaning of the parameters. However 
the Initial values of the parameters should still preferably have the order of magnitude predicted by the formulas above 
so that the parametric model optimization is more likely to determine an absolute minimum and not a local minimum. 
[0064] This transmission line model for twisted pair is based on electromagnetic theory, so it is physical, white-box 
^ model (as opposed to black box), that contains as much a priori knowledge as possible. The model Includes the skin 
effect and proximity effect (up to 2 terms of the series solution) and is considered to be suitable for the best mode of 
operation of the present invention for twisted pair cables. The present Invention is not however, limited to this method. 

B. Scattering Parameter and Transfer Function H of the Access Network 

35 

[0065] In the following the use of the basic parametric model will be described with reference to a lihiited number of 
useful network topologies. 



1) Single Line 

[0066] Figure 1 depicts an equivalent electrical circuit of an access network consisting of a single line. Initial values 
of five parameters, a^ ag a^ a^ %, of the basic parametric model need to be determined and then optimized. The aim 
of the modelling exercise is to obtain an expression for ( 1 -port model, used for the parameter estimation) and H, 
the transfer function ( 2 -port model, needed for capacity estimation and validation). 
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Pa = 
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(13) 



P/ = 



^1 + ^ 



(14) 
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vf = 






with / = 1 ,2 
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[0067] Treating the transmission line as a single port results in: 

[0068] Z, = oo, = 1 , (assumption of an open -circuited end of the transmission line, e.g. at the customer premises 
of a local loop) 



(0069] The following equation defines for a 1 -port model the scattering parameter S,,: 




■2yl 

s - 

^11 “ — 57 / 

1 -pge 



(17) 




(18) 



Z- = 0, p^ = -1.T =2 



H = 






,T/ 



1+P/e 



-2y/ 



(19) 



[0070] T rvc I'^osfo* funaon under consideration starts at the beginning of the line and ends at the load (no generator), 
the iransmt PSD (power spectral density) mask is defined at the location of VI in Fig. 1 , not at the location of Vg. 

2} Homogonoous L cop with 1 Bridged Tap 

[0071] Figure ? depicts an equivalent electrical circuit of an access network consisting of a loop with a bridge tap. 
Initial values of swen parameters, 32 % % ^ % ot the parametric model need to be determined and then 

oplimi/cd The am of the modeling exercise is to obtain an expression for the scattering function and H, the 
transfer function 




( 20 ) 




( 21 ) 



Zf=oo,Pi=1, (assumption open-circuit at the end) 
Tho following oquation defines for a 1-port model the scattering parameter 



5m 
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= 0 . Pp = 1 = 2 (assumption PSD mask defined at V^) 



2-port model: transfer function H 

5 



10 



// = 






+ p,e 



2y. 



H- e 



-2yj 



+ +£-2n/.w,) _p^^-2y(/,+^,) ^3^^^~2r(i,+/,^2/,) 



( 23 ) 



3. Cascade of 2 Line Sections 



[0072] Figure 3 depicts an equivalent electrical circuit of an access network consisting of a cascade of two lines. 
Initial values of ten parameters a^^ a^^ a ^2 ^ ^32 ^42 %2 parametric model need to be determined 

15 and then optimized. The aim of the modeling exercise is to obtain an expression for the scattering function and 
H, the transfer function. 



2/ = «>, P/ = 1 , (assumption open-circuit at the end) 
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P5 ~ « x5 ~ ^ 
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(Yi ^ +Y2^) 



j '^Yi^l “^Y2^2 ~2 (Yi ^1+Y2(2) 

1+p2© +Pl^ 



(25) 



II. ESTIMATORS 

35 

A. Cost function • _ . , 

[0073] In the following one cost function will be described for use in optimizing the parameters of the parametric 
model, but the present invention is not limited thereto. It is called maximum likelihood for the output error model. This 
4c type of cost function has useful properties such as the fact that estimated parameters are not biased, 



45 
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( 26 ) 



M = number of frequency samples 
N = number of parameters 

50 [0074] This first cost function in S,, is important for the general method in accordance with the present invention 

The following second cost function in H is only used for validation puiposes. 
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B. Minimizer 



[0075] The optimization of parameters is carried out by minimizing a difference represented by the cost function. 
The preferred minimizer is from Levenbergh-Marquardt but the present invention Is not limited to this type of optimization 
5 algorithm, nor to any parrticular values for the starting values. 

[0076] As an example, the following starting values may be selected: 

a = 0.25x10*^ mm 
0= 1 *'10*3 mm 

^0 0=5.8 • 1 0^ S/m assumed for copper conductors, 

p = uq for PE insulated cables 
= 2.26 £q for PE insulated cables 
text = mu/(2*pi)*log(D^2/radius^2); 

C = 1/Lext*epsilon*mu; 

15 

The minimum of the cost function as a function of the model parameters may be found by means of an iterative opti- 
myation algorithm. There is a large variety of optimization routines available which are known to the skilled person. 
The Levenborg-Marquardt (LM) method is one of the preferred ones. It combines the Gauss-Newton and gradient 
descent procedures. The Gauss-Newton method is fast, because it makes use of second-order derivatives, but it does 
^ not always converge. This problem of convergence is solved in the Levenberg-Marquardt algorithm by addition of the 
gradient descent method. The Levenberg-Marquardt algorithm only searches for a local minimum. In order to find the 
global minimum, the initial values of the model parameters have to be close enough to the exact values, which corre- 
spond to the global minimurn of the cost function. In accordance with an aspect of the present invention achieving he 
absolute m.nimum of the cost function is assisted by selecting a parametric model which Is at least quasi-realistic so 
^ that ro liable initial values for at least some of the coefficients can be estimated. Preferably, the parametric m*odel should 

bo chosen such that a majority of the coefficients (more than 50%, more preferably 60% or more and most preferably 
70°o or more) can be estimated based on physical characteristics of the network to be qualified. 

Ill MEASUREMENTS 
30 

A. Measurement Setup 

[0077] These are the settings used for the network analyzer measurements for determining (for the estimation) 
and H (for the validation). 

^ [0078] In the measurement step the reflectogram of the network is measured for a specific excitation signal. A Hewlett 

Network Arialyzer type HP4195A was used with a start frequency of 3*4.3125=12.9375 kHz (tone 3) and a 
stop frequency of 1 .104000 MHz (tone 256). The RBW was 100 Hz, sweep time 2 min, and 401 points were recorded. 
Four measurements were performed of and H in each case to be able compute mean and variance values. These 
tatter values are used as weighting factors in the MLE cost function (see equation 26) 

[0079] The transmission lines used were cables, e.g. as supplied by Belgacom (Belgium), having a conductor diam- 
eter of 0.5 mm, a PE insulation, and 20 pairs in the form of quads; as supplied by BT (UK), having a conductor size of 
0.5 mm, insulation of PE, 20 pairs; as supplied by FT (France), having conductor size of 0,4 mm, PE Insulation, 8 pairs, 
quads, the North Hill balance was 50 Ohm to 135 Ohm 

B, Resutts 



[0080] Only the magnitude is shown, although the phase has been estimated as well. The magnitude is however 
most important for the capacity estimation 

[0081] The results are shown for 3 topologies, but the patent is not limited to these 3 topologies. 

LEGEND for graphs in Figs. 4, 6, 8 showing S^i on the y axis is : 

• (dotted) : measured values as shown in upper portion of figures 
- (full line): estimation as shown in upper portion of figures 
-X- (full line, x); error as shown in the middle of the figures 

-4- (full line for Figures 6 and 8 only, +): 95 % confidence bound as shown in lower protion of the figures 

Remark: in the case that there are only stochastic errors present, the error should be below the 95 % confidence bound. 
Where this not the case this means that there are model errors present- 
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LEGEND for graphs in Fig. 5, 7, 9 having H as the y axis is: 

• (dotted) ; measured values as shown In upper portion of figures 

- (full line): estimation via (generally so close to the dots that It cannot be distinguished from these at this scale) 

5 - (dashed line in Figures 7 and 9): estimation via H (generally so close to the dots that it cannot be distinguished 

from these at this scale) 

-X- (full line, x): error in the estimation via in the middle of the figures 

- o- (dashed line in Figure 5, full line in Figures 7 and 9, o): error in the estimation via H in the middle of the figures 

(full line only in Figures 7 and 9, +): 95 % confidence bound in the lower portion of the figures 

10 

1) Single Line 

[0082] In Fig. 4 the results for the scattering factor both model and measurement and in Fig. 5 for the transfer 
function H for a single line, both measurement and model are shown. 

15 

2) Loop with a bridged tap 

[0083] In Fig. 6 the results for the scattering factor S.,i, both model and measurement and in Fig. 7 for the transfer 
function H, both measurement and model are shown. 

[0084] All segments of cable were BT cable with conductor size of 0.5 mm, first segment was 900 m, second segment 
700 m and bridged tap 200 m. 

3) Cascade of 2 tines 



[0085] In Fig. 8 the results for the scattering factor , both model and measurement and in Fig. 9 for the transfer 
function H, both measurement and model are shown. 

[0086] A BT cable of 500 m and conductor size 0.5 mm was cascaded with an FT cable of 400m with a conductor 
siz of 0.4 mm. 

For all the 3 topologies it is clear that the 3 transfer functions (measured, estimated via 1*port measurement 
estimated via a 2-port measurement H) match very well. 

IV Capacity estimation 

[0087] The capacity estimation is done using the formula based on Shannon's theory given as equation 1 above and 
the estimated transfer function H. The noise PSD needs to be determined, e.g. by measurement. At one end of the 
transmission line, e.g. at the CO for an access network, the noise PSD (for the upstream capacity) can be measured. 
At the other end of the line, e.g. the customer premises in an access network, the noise PSD can be estimated for 
example by use of a model for the noise PSD at this end (for the downstream capacity). 

V Implementation 



[0088] The present invention may be implemented on a computing device e.g. a personal computer or a work station 
which has an input device for loading the details of the network whose capacity is to be estimated, e.g. its topology as 
well as any other parameters required by the various models and methods defined above in the description as well as 
in the attached claims. The computing device may be in the form of a card mounted microprocessor which may plugged 
into a network card slot of a telecommunications network element such as the 7300 Alcatel DSL Subscriber Access 
Multiplexer supplied by Alcatel NV, Antwerp Belgium . It could also be integrated in a handheld testing device used by 
a technician in the field or by testing equipment next to the Subscriber Access Multiplexer. The computing device is 
adapted to run software which carries out any of the methods in accordance v/ith the present invention. The computer 
may alternatively be a server which is connected to a data communications transmission means such as the Internet, 
a Local Area Network or a Wide Area Network. A script file including, for example, the details of the topology of the 
network and/or physical constants relating to the materials used in the construction of the transmission line(s), may 
be sent from one near location, e.g. terminal, to a remote, i.e. second location, at which the server resides. In addition 
details of the reflectogram for the relevant network may be sent to the server from the same near location or from 
another location. The server receives all this data and carries out a method in accordance with the present invention 
and outputs back along the communications line useful data to a near terminal, e.g. one or more optimized coefficients 
of a transmission line model, a transfer function of the network, a capacity figure, or a yes no answer as to whether 
the capacity of the network is greater or less than a predetermined target figure. 
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[0089] Fig. 10 is a schematic representation of a computing system which can be utilized in accordance with the 
methods and systems of the present invention. A computer 1 0 is depicted which may include a video display terminal 
14, a data input means such as a keyboard 16. and a graphic user interface indicating means or pointer such as a 
mouse 18. Computer 10 may be implemented as a general purpose computer. 

[0090] Computer 10 Includes a Central Processing Unit (“CPU") 15, such as a conventional microprocessor of which 
a Pentium III processor supplied by Intel Corp. USA is only an example, and a number of other units interconnected 
via system bus 22. The computer 1 0 includes at least one memory. Memory may include any of a variety of data storage 
devices known to the skilled person such as random-access memory ("RAM”), read-only memory ("ROM”), non-volatile 
read/write memory such as a hard disc as known to the skilled person. For example, computer 10 may further include 
random-access memory ("RAM”) 24, read-only memory ("ROM") 26, as well as an optional display adapter 27 for 
connecting system bus 22 to an optional video display tenrilnal 14, and an optional input/output (I/O) adapter 29 for 
connecting peripheral devices (e.g., disk and tape drives 23 such as for example a CD-ROM reader) to system bus 
22. Video display terminal 14 can be the visual output of computer 10, which can be any suitable display device such 
as a CRT-based video display well-known In the art of computer hardware. However, with a portable or notebook- 
based computer, video display terminal 1 4 can be replaced with a LCD-based or a gas plasma-based flat-panel display. 
Computer 1 0 further includes user interface adapter 30 for connecting a keyboard 1 6, mouse 1 8, optional speaker 36, 
as well as allowing optional physical value Inputs from physical value capture devices 40 of an external system 20. 
The devices 40 may be any suitable equipment for capturing physical parameters of the network or parameters of the 
various models required in the execution of the present invention. These capture devices may also include a stimulus 
and a measurement device for inputting a test wave and for measuring the response of a network or parts thereof, e. 
g. a reflectogram. Additional or alternative devices 41 for capturing physical parameters of an additional or alternative 
external system 21 may also connected to bus 22 via a communication adapter 39 connecting computer 10 to a data 
network such as the Internet, an Intranet a Local or Wide Area network (LAN or WAN) or a CAN. The term “physical 
value capture device" can also includes devices which provide values of parameters of a network or networks, e.g. 
topologies, or for instance a library of candidate networks or candidate network topologies. 

[0091] Computer 1 0 also includes a graphical user interface that resides within a machine-readable media to direct 
the operation of computer 1 0. Any suitable machine-readable media may retain the graphical user interface, such as 
a random access memory (RAM) 24, a read-only memory (ROM) 26, a magnetic diskette, magnetic tape, or optical 
disk (the last three being located in disk and tape drives 23). Any suitable operating system and associated graphical 
user interface (e.g. Microsoft Windows) may direct CPU 1 5. In addition, computer 1 0 includes a control progreun 51 
which resides within computer memory storage 52. Control program 51 contains instructions that when executed on 
CPU 15 carries out the operations described with respect to the methods of the present invention. The instructions 
may be obtained by writing a computer program in a suitable language such as C or C-m- for execution of any of the 
methods in accordance with the present invention and then compiling the program so that it executes on a computing 
device. 

[0092] Those skilled In the art wijl appreciate that the hardware represented in_FIG. 10 may vary for specific appli- 
cations. For example, other peripheral devices such as optical disk media, audio adapters, or chip programming de- 
vices, such as PAL or EPROM programming devices well-known In the art of computer hardware, and the like may be 
utilized in addition to or in place of the hardware already described. 

[0093] In the example depicted in Fig. 1 0, a computer program product (i.e. control program 51 for executing methods 
in accordance with the present invention comprising instruction means in accordance with the present invention) can 
reside in computer storage 52. The instructions (e.g., computer readable code segments in storage 52) may be read 
from storage into RAM 24. Execution of sequences of instructions contained in the RAM 24 causes CPU 15 to perform 
the process steps described herein. In alternative embodiments, hard-wired circuitry may be used in place of or in 
combination with software instructions to implement the invention. Thus, embodiments of the invention are not limited 
to any specific combination of hardware circuitry and software. Accordingly, the present invention may take the form 
of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and 
hardware aspects. 

[0094] Furthermore, the present invention may take the form of a data carrier medium (e.g. a computer program 
product on a computer-readable storage medium) carrying computer-readable program code segments embodied in 
the medium. The terms "carrier medium" and "computer-readable medium" as used herein refer to any medium that 
participates in providing instructions to a processor such as CPU 15 for execution. Such a medium may take many 
forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non-volatile media 
include, for example, optical or magnetic disks, such as a CD-ROM or a storage device which is part of mass storage. 
Volatile media includes dynamic memory such as RAM 24. Transmission media include coaxial cables, copp r wire 
and fiber optics, including the wires that comprise a bus within a computer, such as bus 22. Transmission media can 
also take the form of acoustic or light waves, such as those generated during radio wave and infra-red data commu- 
nications. 
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[0095] Common forms of computer-readable media include, for example a floppy disk, a flexible disk, a hard disk, 
magnetic tape, or any other magnetic medium, a CD-ROM, any other optical medium, punch cards, paper tapes, any 
other physical medium with patterns of holes, a RAM, a PROM, an EPROM, a FLASH-EPROM, any other memory 
chip or cartridge, a carrier wave as described hereafter, or any other medium from which a computer can read. 
[0096] These various forms of computer readable media may be involved in carrying one or more sequences of one 
or more instructions to processor 15 for execution. For example, the instructions may initially be carried on a magnetic 
disk of a remote computer. The remote computer can load the instructions into its dynamic memory and send the 
instructions over a telephone line using a modem. A modem local to the computer system can receive the data on the 
telephone line and use an infrared transmitter to convert the data to an infra-red signal. An infra-red detector coupled 
to a bus can receive the data carried in the infra-red signal and place the data on the bus. The bus carries data to main 
memory, from which a processor retrieves and executes the instructions. The instructions received by main memory 
may optionally be stored on a storage device either before or after execution by a processor. The instructions can also 
be transmitted via a carrier wave in a network, such as a LAN, a WAN, or the Internet. However, it is Important that 
those skilled In the art will appreciate that the methods of the present invention are capable of being distributed as a 
program product in a variety of fornis, and that the present Invention applies equally regardless of the particular type 
of signal bearing media used to actually carry out the distribution. Examples of computer readable signal bearing media 
include: recordable type media such as floppy disks and CD ROMs arid transmission type media such as digital and 
analogue communication links which may be used for downloading the computer program product. 

VI Method 

[0097] A method in accordance with the present invention will be described with reference to a personal computer. 
A computer program written to cany out the methods of the present invention is started in the usual way. For instance, 
a pop-up message appears providing an input form for specifying the network topology to be used. Alternatively or 
additionally the possibility for using predefined topologies may be provided, e.g. by allowing selection from a drop- 
down menu. Alternatively or addrtionalty an option may be provided when the topology is not known. In this case the 
program will make use of a set of predefined typical topologies and will attempt to fit a measured reflectogram to each 
of the models. The best fit is assumed to represent a topology which is close to the actual topology and is used for the 
capacity estimation. 

[0098] After having entered the necessary data, for example network topology, insulating and conductor materials 
of a twisted pair cable, conductor spacing of the cable, etc. the reflectogram for the network may be entered. This may 
be done by loading from a peripheral device, e.g. from a CD-ROM read by a CD-ROM reader, or by downloading it 
from mass storage, e.g. from a server located on a LAN or via the Internet. Alternatively, the computer may be connected 
to a suitable reflectogram capture device, e.g. via a USB interface, for direct measurement and loading of the reflec- 
togram. The reflectogram may be obtained from a single location in the network, e.g. from the connection point of a 
local loop to a central office. The reflectogram will generally be in analogue form and is converted into digital signals 
by an analogue/digital converter. The program then optimizes the coefficients of the parametric model, e.g. by mini- 
mizing a cost function relating to a difference between the measured reflectogram and a computed transfer function 
for the model. Using the optimized parametric model a transfer function of the network is determined. From this transfer 
function a capacity for the network can be estimated. 

[0099] The output of the computer program can be any or all of the following: 

• a set of optimized parameters and/or a network topology which provides the best fit 

• a transfer function of the network 

• an estimated reflectogram based on the optimized parametric model 

• a comparison between the optimized reflectogram and the measured reflectogram 

• an estimated capacity for the network 

• a yes/no answer or a probabiity that the network will support a capacity above or below a certain target capacity 

• a message that the reflectogram cannot be modeled with a reasonable accuracy (due to unusual topology, unusual 
materials) 

• a length of the network which can support a certain capacity or a variety of lengths for a variety of capacities. 

list of symbols 

C = channel capacity 
P = received power spectral density 

S = transmitted power spectral density OR scattering matrix 
N = noise power spectral density 
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H = transfer function 
7 = propagation constant 
2q = characteristic impedance 
/= length 

5 G = electrical conductivity OR standard deviation 

p = magnetic permeability 
e = electrical permittivity 
D - distance between conductor centers 
a = radius conductors 
10 = physical model parameter 

s = Laplace frequency 
<D = angular frequency 
Zg = per-unit-length serial impedance 
Yp = per*uniMength parallel admittance 
15 y = auxiliary function 

Jf= Bessel function 
= voltage 
/y = current 

= forward (+) and backward (-) travelling voltage wave 
^efi- reference impedance 
py = reflection coefficient 
^ transmission coefficient 
Sfj = scattering parameters 
T = transmission (ABCD) matrix 
^5 A, B, C, D= elements of ABCD matrix 

Z,„ = input impedance 



Claims 
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1 . A method of modeling a netwofic, comprising the steps of: 

carrying out a 1 port measurement on the network by inputting an excitation signal at one port of the network 

and recording the results at the port; 

generating a 1 port parametric model of the network; 

, . optimising the values of the parameters of the 1 port parametric_rnodeLby reducing , the difference. between 

the results of the measurement step and results calculated using the 1 -port parametric model and the excitation 
signal. 

40 2. The method according to claim 1 , wherein the reducing step is earned out by minimizing a cost function by adapting 

the values of the parameters of the 1 port parametric model, the cost function representing a quantification of the 
difference. 

3. The method according to claim 1 or 2, further comprising the steps of: 

45 

generating a 2-port parametric model using the optimized parameter values to provide a 2-port model of the 
network. 

4. The method according to claim 3, further comprising the step of limiting the 2-port model to the transfer function. 

50 

5. The method according to claim 4 further comprising the steps of: 

measuring a first noise PSD at one port of the network and estimating or measuring a second PSD at the other 
port; and 

using the transfer function and the first and second noise PSD to estimate the capacity of the network. 

6. A system for modeling a network, comprising: 
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means for carrying out a 1 port measurement on the network by inputting an excitation signal at one port of 

the network and recording the results at the port; 

means for generating a 1 port parametric model of the network; 

means for optimizing the values of the parameters of the 1 port parametric model by reducing the difference 
between the results of the measurement step and results calculated using the 1-por1 parametric model and 
the excitation signal 

7. The system according to claim 6, wherein the means for reducing comprises means for minimizing a cost function 
by adapting the values of the parameters of the 1 port parametric model, the cost function representing a guanti* 
ficatlon of the difference. 

8. The system according to claim 6 or 7, further comprising : 

means for generating a 2-port parametric model using the optimized parameter values to provide a 2-port 
model of the network. 



9. The system according to claim 8, further comprising: means for limiting the 2-port model to the transfer function. 

10. The system according to claim 9 further comprising: 

means for measuring a first noise PSD at one port of the network and estimating or measuring a second PSD 
at the other port; and 

means for using the transfer function and the first and second noise PSD to estimate the capacity of the 
network. 



11 . A data carrier medium carrying one or more computer readable code segments for controlling a processing device 
to carry out a method in accordance with any of the claims 1 to 5. 

12. A computer program product for execution of any of the methods in accordance with any of the claims 1 to 5 on 
a computer system. 

13. A method for modeling a representation of networ1<, said method comprising: 

transmitting a description of the network from a near location to a remote computing system via a telecom- 
munications network, 

executing on the remote computer system any of the methods of claims 1 to 5, and transmitting to a near* 
location information relating to the response of the network to an excitation signal. 




parametric line model 
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